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Abstract. A thin tube is an n-dimensional space which is very thin in n— 1 di- 
rections, compared to the remaining direction, for example the e-neighborhood 
of a curve or an embedded graph in R" for small e. The Laplacian on thin 
tubes and related operators have been studied in various contexts, with differ- 
ent goals but overlapping techniques. In this survey we explain some of these 
contexts, methods and results, hoping to encourage more interaction between 
the disciplines mentioned in the title. 
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1. Problem 

A thin tube is, roughly speaking, a metric space which is close to a one- 
dimensional space, for example a curve. Often, but not always, it is useful to 
think of a thin tube as a member of a family {ng)o<e<eo where fig is built from 
a one-dimensional space Hq by adding transverse dimensions scaled to size e in a 
prescribed way. Thus, JIq can be thought of as limit (as e — > 0) or model for the 
family {fle)e- 

We will make this precise in various settings. We will consider two kinds of 
tubes: Tubes without ends and tubes with ends. To explain this, we define these 
first in the simple setting of two-dimensional tubes modeled on an interval. See 
Figure 1. 

Tubes without ends: Let / C M be a closed interval and let hi, h2 : I ^ M. be 
continuous functions satisfying hi{x) < h2{x) for all x E I. For e > the tube 
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Figure 1. Tube without ends and tube with ends 

without ends is 

(1.1) f^without ends ^ |(^^^) g . ^ g eh2{x)}. 

Tubes with ends: In addition to the previous data, let X;, Xr be two bounded 
closed sets in with non-empty interior such that Xi is contained in the left half 
plane {x < 0} and X^ in the right half plane {x > 0}. A tube with ends is obtained 
by attaching Xi and X,., scaled by a factor e, to the left and right of a tube without 
ends. Thus, for / = [a, 6] 

(1.2) nf*^ """^^ = f7f ^^'^ U [(a, 0) + sXi] U [{b, 0) + eXr]. 

where a + eX := {(a + ex,ey) : {x,y) E X}. We assume Xi,Xr are such that 
j-^with ends jg connected. Note that the ends are scaled by e in all directions while 
the tube is only scaled in the y-dircction. This will make a substantial difference 
in the analysis. 

We simply write f2e if we do not want to specify which kind of tube we consider. 



The class of problems we will discuss concern the Laplace operator, A = 
— + on rig. We will be mostly interested in the spectrum, that is, the 
solutions of 

(1.3) Au = Xu, u:n.~^C,XeR 

which in addition satisfy boundary conditions, for example Dirichlet {u^dn^ — 0) 
or Neumann {d„u\QQ^ = where dn is the derivative normal to the boundary) 
conditions. We assume boundedness of fie and sufficient regularity of its boundary 
to ensure that the spectrum is discrete with finite multiplicities. Therefore, for each 
£ > we obtain a sequence 

Ai(e) < A2(e) < > oo 

of eigenvalues, where each eigenvalue appears as often as the dimension of its 
eigenspace indicates. Corresponding eigenfunctions will generally be denoted Uk^e- 

Main Problem. Describe the asymptotic behavior of the eigenvalues Xk{£), of 
spectral invariants constructed from these, and of the eigenspaces as e — > 0. 

We now describe some other sorts of thin tubes, which will be referred to 
throughout the article. 

• Manifolds and higher dimensions: Instead of plane domains one 
may consider domains in R" and more generally Riemannian manifolds. 
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with the Laplace-Bcltrami operator.^ Thus, for some n > 2 we start 
with an (n — 1) dimensional manifold Y and two n-dimensional manifolds 
Xi,Xr. All manifolds are compact and Riemannian and may have (piece- 
wise smooth) boundary.^ Also, we are given isometrics ji, jr from Y to 
subsets of the boundary of Xi, Xr- For e > and a Riemannian manifold 
{X,g) we denote by eX the Riemannian manifold {X,e^g), that is, all 
lengths are multiplied by e. For an interval I = [a,b] we then set 



{the thin tube manifold with cross section Y and ends Xi,Xr) where the 
union is disjoint and the quotient means identiiying a x Y with ji{Y) C 
dXi using ji and bxY with jr(^) C dXr using jV- 

Note that if Y is an interval and Xi,Xr C M? then this reduces to 
(1.2) with hi and /12 constant. It is straightforward to give a manifold 
setting where the cross-section has varying thickness. 

Note also that this setting includes the e-neighborhood of a straight- 
line segment in M" (that is, points of distance at most e from the line), 
as well as the boundary of this neighborhood (sometimes referred to as 
sleeve). In the former case, Xi,Xr are half balls and Y is the (n — 1)- 
dimensional unit ball, and in the latter case Xi,Xr are half spheres and 
F is a sphere. 

The Laplace(-Beltrami) operator on a space X, with boundary con- 
ditions given in the context, will be denoted by Ax- We always take it 
positive, i.e. in M" it is 



• Curved and twisted tubes: Instead of modeling fig on a straight 
line, one may consider domains obtained as £-neighborhood of an arbi- 
trary simple curve. More generally, instead of a circular cross-section one 
may look at arbitrary cross-sections. The cross section may also rotate 
when progressing along the curve. Such tubes may be considered with or 
without ends. 

• Tubes modeled on graphs: An important generalization of tubes 
with ends is obtained by replacing the single line segment by a metric 
graph. Here, a graph is given by sets V of vertices and E of edges, where 
each edge connects two vertices (or one vertex with itself). Write e ^ v 
if edge e has vertex v as an endpoint. A metric graph is a graph together 
with a fimction I : E ^ (0, 00). /(e) is to be thought of as the length 
of e. Suppose we are given an (n — l)-dimensional manifold Fe for each 
edge e and an n-dimensional manifold Xy for each vertex v, as well as 
isometrics (gluing maps) of Y,, with a subset of the boundary of Xy for all 
e ^ V, without overlaps (i.e., the images of the gluing maps corresponding 



Instead of the Laplace-Bcltrami operator, one may also consider other operators relating to 
the tube geometry, for example the Dirac operator. 

■^for an exact definition of this, see for example Hassell and Zelditch (2004); for our purposes 
one may assume that the data are such that the spaces defined below are smooth manifolds with 
boundary. 



(1.4) 



= [eXi U (/ X eY) U eX,] 
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to a single vertex are disjoint). Then fl^ is defined by gluing cylinders of 
length Z(e) and cross section eYe to eX^ for all pairs e ~ w. 

A special case is the e-neighborhood of a graph embedded in R" with 
straight edges, and the boundary of this neighborhood. Of course, one 
may also consider curved and twisted edges. 

The additional combinatorial structure of the graph suggests a refine- 
ment of the Main Problem in this context: 

Quantum Graph Problem. How is the limiting behavior of the 
spectral data on fig related to spectral data of some operator on the metric 
graph itself? 

This question naturally leads to quantum graphs, to be introduced 
below. 

• Noncompact tubes: So far, all tubes were compact. Instead, one may 

consider tubes modeled on the real line or half line, or more generally 
on graphs having some edges of infinite length (which then have only one 
endpoint). Such spaces arise in their own right and also as certain rescaled 
limit objects of families of compact tubes, see equation (3.21). Apart from 
this, however, we focus mostly on compact tubes in this article. 
The goal of this article is to explain some common ideas and methods used in 
different areas where thin tubes occur. We also explain a very limited selection of 
mostly recent results, mainly from mathematical physics and spectral geometry. 

Thin tubes as defined above are a special class of singular degenerations. There 
are many other kinds of singular degenerations. For example, one may keep the 
vertex neighborhoods fixed while only shrinking the tubes between them (this is 
sometimes called a handlebody), or one may remove a small ball inside some fixed 
domain. There is a large body of literature on such problems. We just mention the 
book by Maz'ja, Nazarov, and Plamenevskij (1991). See also the discussion in 4.4. 

In addition to the references given below, we would like to mention the collection 
by Bcrkolaiko ct al. (2006) which has many articles on related questions. 

The author acknowledges the hospitality and support of the Isaac Newton In- 
stitute in Cambridge during the program 'Analysis on Graphs and its Applications' 
in 2007, where part of this work originated. 

1.1. Rescaling to long cylinders. For a proper understanding of many tube 
problems it is important to consider tubes rescaled to fixed width. Denote 

(1.5) N = e-\ X^=e-^ne. 

That is, all lengths in fig are stretched by the factor e~^. This gives a tube of fixed 
width and with fixed ends (independent of e) and length on the order of = 
(assuming fig is given as in (1.2) or (1.4), and similarly for tubes modeled on a 
graph). The eigenvalues and eigenfunctions of the Laplacian on are^ 

(1-6) Xk{N) = e'^Xk{e), Uk,N{x,y) = Uk,e{£x,ey). 

Note that the limit e — > corresponds to A'' — > oo, and this yields a limit object 

X°° different from ilo- While ilo is just the one-dimensional skeleton of the tube, 
X°° has the same dimension as the tube and also knows about the shape of the 

^We will denote rescaled quantities and variables with a bar, for example A, u for eigenvalue 
and eigenfunction, and x,y for coordinates on X". 
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cross sections and of the ends. In this sense, X°° retains more information about 
the tubes than Oq. On the other hand, in the case of graphs, fJo knows about the 
graph (connectivity) structure while does not. The precise definition of X°° 
will be given in (3.21). 

In fact, in part of the literature it is common (and more natural) to state 
problems and results in the 'long cylinder' setting, especially in the global analysis 
literature. 

In this article I will state most results in terms of thin tubes, so that they can 
be compared more easily. 

2. Origins 

2.1. Mathematical Physics. There arc a number of physical systems which 
are quasi one-dimensional in the sense that their dimensions in one direction are 
much larger than their dimensions in transversal directions, or more generally that 
several such onc-dimcnsional subsystems arc coupled in a graph-like structure. The 
quantum mechanics of such systems is then most easily modeled by the underly- 
ing graph (i.e. neglecting all transversal directions) with a Schrodinger operator on 
the edges coupled by transition (boundary) conditions at the vertices (a so-called 
quantum graph). A better, still simplifying model, is to consider thin tubes around 
the graph edges and a Schrodinger operator, for example the Laplacian on such 
tubes. Examples of such systems are: carbohydrate molecules, where the graph 
structure is given by the nuclei and their bonds - which are assumed fixed -, and 
the goal is to analyze the behavior of the electrons; this in fact seems to be the 
earliest occurrence of such a model, see Rucdenbcrg and Scherr (1953). Another, 
more recent kind of such quasi one-dimensional systems that need to be analyzed 
quantum-mechanically are so-called nanotubes or quantum wires, or highly inte- 
grated circuits. See the nice survey by Kuchmcnt (2002) for more background and 
many references, and the classical article by Exner and Seba (1989b) which sparked 
a lot of activity in this field. 

Confining a quantum-mechanical particle to a small neighborhood of a graph- 
like structure may be achieved also by other means than by a thin tube: For 
example, one may consider the Schrodinger operator on M" with a potential a'^V 
added, where a is a large parameter and V is non-negative, and vanishes precisely 
on the graph. This leads to closely related questions and results. See for example 
Smilansky and Solomyak (2006), where a transversal harmonic oscillator is added 
to the quantum graph operator, and Ftoese and Herbst (2001) (in the case of 
submanifolds instead of graphs). 

2.2. Global analysis. Global analysis studies how analytic objects associated 
with certain partial differential operators on a manifold (e.g. spectrum, dimension 

of kernel) arc related with global, e.g. topological, properties of the manifold. This 
includes for instance Hodge theory (where the cohomology of a compact manifold 
is calculated in terms of the kernel of the Laplace-Beltrami operator on differential 
forms) and index theory (where the difference of the dimensions of the solution 
spaces of a system of PDEs on a compact manifold and of its adjoint is computed 
by topological means). A typical procedure is this: one starts with a smooth 
compact manifold X, then chooses a Ricmannian metric g on X, this defines the 
Laplace-Beltrami operator Ax (or, in the presence of additional structure, Dirac 
operators), then one forms certain expressions from the spectrum of Ax, and finally 
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these are shown to be independent of the choice of metric, and hence (diffcrential- 
)topological invariants. This happens in Hodge theory, and another instance of this 
is the analytic torsion, which is a combination of determinants of the Laplacians 
(introduced below in 3.3.2) on A:-forms, for A; = 0, . . . ,dimX. The determinant of 
the Laplacian itself is only a spectral invariant, i.e. it depends on the eigenvalues 
(and hence on the metric). There are also other spectral invariants (for example, 
the ry-invariant for Dirac-operators) which appear in global analysis contexts. 

It is then a basic problem to understand how these topological invariants change 
under modifications of the manifold: For example, if one cuts the manifold into 
two pieces by removing a hypersurface Y then one would like to know how the 
invariant of the whole manifold relates to the invariants of the pieces (one seeks a 
gluing formula). 

The same question can be asked for spectral invariants. This is interesting in 
its own right, but since the topological invariants are defined in terms of these, an 
answer to the spectral question will also answer the topological question. 

Surprisingly, a useful technique for obtaining such gluing formulas is to consider 
a family of Riemannian manifolds X'^ obtained by gluing a cylinder'' [—N, N]xY 
into X in place of y, and studying the asymptotic behavior of the spectral invariants 
of X^ as A?^ — > 00, in terms of those of Y and of the ends Xi, Xr- An example 
of this will be given in 3.3.2. As manifold, X^ is the same as (more precisely, 
diffeomorphic to) X. It is only the metric which is different. Since the topological 
invariants are independent of the metric, hence of A'', this asymptotic behavior 
determines the behavior also for A'' = 0, thus providing the desired gluing formula. 

2.3. Spectral geometry. Spectral geometry studios how the eigenvalues and 
eigcnfunctions of the Laplacian on a space depend on the geometry of the space, 
and which properties they have independent of the geometry. In this context, thin 
tubes sometimes serve as a special class of spaces for which questions that are very 
hard in general may be answered by asymptotics methods. We give two examples 
of this: 

The nodal line conjecture 

W; consider the Dirichlet Laplacian on a bounded plane domain O. It is well-known 
that the first (smallest) eigenvalue is simple, that the corresponding eigenfunction 
does not change sign in O and that consequently any eigenfunction u for the second 
smallest eigenvalue must have a non-empty nodal lineAf := u~^{0). Also Courant's 
nodal domain theorem implies that TV divides f2 into two connected components, 
and this implies that it is a smooth curve. Payne (1973) conjectured that Af must 
touch the boundary of il, that is, it cannot be a closed curve in the interior of fl.^ 
Jerison (1991) approached the problem in the case of convex domains ft as 
follows: By scale invariance, one may assume that the diameter of O (the largest 
distance between two points in Q) is one. Let s be the inradius of fl, that is the 



to obtain a smooth metric, one has to assume that the metric on X is of product type, i.e. 
dx'^ + ffy , in a neighborhood of Y 

®The state of the nodal line conjecture is as follows; Melas (1992) proved it for convex 
plane domains. Hoffmann- Ostcnhof, Hoffmann- Ostcnhof, and Nadirasfivili (1997) showed that it 
cannot hold if one allows the domain to be multiply connected. Jerison (1995b) proved it for 
convex domains in R" with diameter one and small inradius. An interesting higher dimensional 
result was also obtained by Freitas and Krejciffk (2007). An open problem is whether the nodal 
line conjecture is true for plane simply connected domains. 
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maximum radius of a ball contained in 17. It is not hard to show that f2 is of the 
form (1.1), where hi, h2 are uniformly boimded with max(/i.2 — hi) of order one. 

Therefore, if the inradius is small then this is a case of thin tubes, and therefore 
some of the techniques described below are applicable in principle. However, a 
major complication arises from the fact that the first and second derivatives of hi , /12 
cannot be controlled beyond what follows from convexity. In spite of this, Jerison 
manages to prove the nodal line conjecture for convex domains with sufficiently 
small inradius. See 3.2.2 for more on this. 

Multiplicity of the second eigenvalue 

A basic question in spectral geometry is: Given a compact manifold M without 
boundary and an increasing sequence of non-negative real numbers, is there a Rie- 
mannian metric on M such that the given sequence is precisely the eigenvalue se- 
quence of the Laplace-Beltrami operator? This is probably intractable, so one may 
simplify the question by just prescribing a finite number, say the first N eigenvalues 
(always counted with multiplicities). 

Assuming M is connected and without boimdary, a first condition on the se- 
quence is that the first eigenvalue be zero and have multiplicity one. In two di- 
mensions there is another well known obstruction: The multiplicity of the higher 
eigenvalues is bounded in terms of the Eulcr characteristic of M, see Cheng (1976). 

In dimensions n > 3, Colin dc Vcrdicrc (1986) showed that there is no such 
obstruction for the multiplicities of the second eigenvalue^. The idea of the proof 
is as follows: Given N E N. embed a complete graph with N vertices in M (here 
one needs that dimAf > 3, if is large). On a neighborhood of this graph in 
M, define the metric to be of the thin tube type H,^ (with all edges of length one 
and disks as cross sections), and extend it outside of this neighborhood so that 
it is small there in a suitable sense. Then the eigenvalues on M will be close to 
those on fig, with Neumann boundary conditions at dfle, and these will be close to 
those of a suitable operator on the graph itself if e is small (in the latter step the 
methods discussed below are used; see below for the definition of the operator on 
the graph). A simple calculation shows that the second eigenvalue of the operator 
on the graph has multiplicity A'' — 1. Hence the second to A^th eigenvalue of M 
must lie very close together. A topological perturbation argument then shows that 
one may change the metric on M to obtain one where the multiplicity of the second 
eigenvalue equals N — 1. 

3. Results 

We will explain a number of results and try to motivate why they should be 
true. The selection is restricted by space as well as the personal taste and the 
limitations of the author. I apologize to anyone (that is, almost everyone) whose 
result is not included. 

We will focus on results directly related to spectral quantities. In some cases, 
there are actually stronger statements in the papers concerning the resolvent, which 
imply the spectral results. See the explanation in 4.2. We emphasize results where 
Dirichlet (or other non-Neumann) boundary conditions are imposed, since this case 



In fact, Colin de Verdiere (1987) shows the stronger result that any sequence of N numbers 
= Ai<A2<-<Ajv can occur. 
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has seen the greatest progress recently. Ideas of proofs of some of the resuhs are 
given in Section 4. 

After introducing the basic principle that guides what one should expect for 
thin tube problems we consider results where the limit object is an interval or a 
curve and then the case of graphs. Clearly intervals are a special case of graphs. 
However, there are many results that are only natural for intervals, or have so 
far only been formulated for these, or are better (for example, sharper) than their 
graph analogues. 

3.1. The basic principle. The basic principle in thin tube problems is that 
for the lower part of the spectrum the first transversal mode dominates the behavior 
of eigenvalues and eigenfunctions. To see this, we first look at the simple example 
of tubes modeled on an interval / and with constant cross section Y. Let 

(3.1) I'D < z^i < • • • the eigenvalues of Ay, with eigenfunctions (po, ifi, . . . , 

with the boundary conditions given at dY. Then the eigenfunctions and eigenvalues 

on = / X eY are^ 

(3.2) Us{x, y) = sm{aix) (Pm(|), A(e) = e~'^Urn + ah = 

for Z e N, m e No (for Dirichlet conditions at the ends of /; replacing sin by cos 
and taking I G No gives Neumann conditions). So for any fc e N, the fcth eigenvalue 
is £"^1^0 + ttfc as soon as s is sufficiently small. This leads us to expect for general 
tubes that 

• for the (rough) behavior of Xkie) for fixed k and £ ^ only the first 
transversal mode (m = 0) should be relevant, and 

• higher transversal modes (m > 0) will be relevant in the behavior of Afc(£) 

if, simultaneously with e ^ 0, one lets fc — > oo on the order of e~^. 

In the case of a plane domain with variable cross section, (1.1), the first eigenvalue of 
the transversal operator at x, —d'^/dy'^ on the interval [h\{x), h2{x)\ with Dirichlet 

conditions, is 7r^/ft(.T)^ with h{T) = h2{x) — hi(x), and correspondingly one expects 
the behavior of low eigenvalues to be dominated by those of the one-dimensional 
Schrodinger operator 

(with boundary conditions at dl corresponding to those of the original problem) 
and the eigenfunctions to be approximately 

7r 

(3.4) Us{x, y) « uo{x) sm——{y - hi{x)), 

en{x) 

for eigenfunctions uq of Pq. 

There are many results that make these expectations precise in one way or 
another. They differ in the generality of the setup (constant or variable width; 
tubes with or without ends; interval, curve or graph as limit object; regularity 
assumptions etc.) as well as in the precision that is aimed at. The desired precision 



'^A remark on notation: If Y is a Riemannian manifold, instead of a domain in R"^-'^, tlien 
y/e does not really make sense. We use it anyway since it is suggestive. The precise meaning is: 
2/ is a point of Y considered with the rescaled metric e^gy, and y/s is the same point in Y, where 
we now use the metric gy. 



THIN TUBES SURVEY 



9 



is motivated by the origin of the problem: In some problems interesting behavior 
occurs already in the leading order asymptotics, while for some one has to go to 
second or third order (in which case it may be natural to consider asymptotics to 
arbitrary order). Having higher order asymptotics may also be useful since in some 
cases it simplifies the proof of the validity of the lower orders. In reverse, this means 
that in irregular problems (where there can be no higher order asymptotics) it is 
technically harder to prove validity of even leading order asymptotics. An example 
of this is (Grieser and Jerison, 1996, 1998), see Section 3.2.2. 

3.2. Intervals and curves as limits: Tubes without ends. As will be 

scion when considering graphs (sec 3.4.3) the c;ase of Neumann conditions is rather 
uninteresting if Q,q is a single curve. Therefore, it is not surprising that all the re- 
sults mentioned here assume Dirichlet boundary conditions. There is one exception 
concerning nodal lines in 3.2.2. 

3.2.1. Domains of variable thickness, I. Friedlander and Solomyak (2007a, b) 
consider a domain (1.1) with hi = 0, and they assume that h = h2 has a unique 
maximum on /, say at the interior point x = 0, and that near a; = it has pih 
order behavior (normalized) 



(3.5) h{x) = 



1-c+xP + OixP+^), x>0, 
l-C-\x\P + 0{\x\P+^), x<0, 



while it is on / \ {0}. Here c± > and p>l. For example, for p = 1 a triangle 
is allowed. One considers Dirichlet boundary conditions at the upper and lower 
boundary. 

What should we expect in this setting? According to the basic principle of 
tube problems one expects the behavior of eigenvalues to be dominated by the 
Schrodinger operator (3.3). This is, in terms of its dependence on the small param- 
eter s, a semi-classical operator (replace e by h), and for these it is standard that 
eigenfunctions concentrate strongly near the miniirnim of the potential V{x) = j^, 
i.e. near a; = 0. So the values of V away from x = should have 'small' in- 
fluence on the behavior of eigenvalues and eigenfunctions. Therefore, we replace 
V{x) =TT^ {\ + 2c±\x\P + 0{\x\P+^) (from (3.5)) on / by tt^ 2tt^c±\x\p on M. The 
first summand gives a contribution of £~^7r^ to the eigenvalue, and the second a 
contribution equal to an eigenvalue of the Schrodinger operator +£~'^G±\x\p , 
where C± = 2tt'^c±. Rescaling a; = e"^, a = shows that these eigenvalues 

equal s~^°'iJ.k where 

(3.6) fik are the eigenvalues of — —-^ + C±\^\p on R. 
Summarizing, one should expect 

(3.7) Afc(£) = £-V+e-%fe + o(£-2«), a = 



p + 2 



This is precisely the result of Friedlander and Solomyak (2007a) (in the case of 
Dirichlet conditions at the left and right ends) and Friedlander and Solomyak 
(2007b) (in the case of Neumann conditions at the left and right ends, and in 
a similar form with the interval / replaced by M). Their main tool is a resolvent 
analysis, see 4.2. They also prove an approximation (3.4) of the eigenfunctions (in 
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L^-norm), whore uq can be replaced by the appropriately rescaled eigenfunctions 
of the operator in (3.6). 

3.2.2. Domains of variable thickness, II: Low regularity. Jerison (1995a) and 
Grieser and Jerison (1996, 1998) consider domains (1.1) but only assume 

• convexity, i.e. h\ is convex and /12 is concave, 

• uniform bounds on /ii, /12; we normalize^ these so 

max-h{x) = 1, h = h2 — h\. 

There are no bounds on the derivatives of h, it is not even assumed to be differen- 
tiable.^ Put differently, one does not fix hi, /12 and then lets £ ^ 0, but rather one 
considers, for any small e, all possible /ii, /12. Note that in this class the distinction 
'with end' - 'without ends' is immaterial since an end (in the convex case, say) 
corresponds to hi/2 having derivative on the order of in an £- neighborhood 
of dl (plus a trivial 1 + Ce rescaling of length). One assumes Dirichlet boundary 
conditions. 

Eigenvalues: Can one expect a two-term asymptotics for the eigenvalues? 
Since in this setting one finds domains of the form (3.5) with arbitrary p, one may 
expect from (3.7) that the answer is no. 

Surprisingly, one may derive from /ii,ft.2 a geometric quantity which almost 
yields a second term for Xk{e)-}^ Jerison (1995a) proves that there are absolute 
constants c, C > such that 

(3.8) ,-V + ^<A,(£)<£-V + ^ 
where I is the biggest number such that 

(3.9) /i > 1 — on an interval C / of length I. 

For example, if / = [0, 1] then ^ = 1, I' = I for a rectangle {hi, h2 constant) while 
£ = £2/3^ J _ jQ^^j £qj, triangle {hi = 0, h2{x) = 1 — x), and it can be shown that 
£ always lies between these values. Compare also with (3.7) where £ w e". i should 
be regarded as a characteristic length scale for Cl^. 

Eigenfunctions: To which precision does the approximation (3.4) hold? More 
generally, to what precision can qualitative properties of be deduced from a 
knowledge of h alone?^^ Normalize by assuming max|Me| = 1, Jerison (1995a) 
shows that 'lives' near the interval I( in (3.9) in the sense that it decays like 
g-d(x)/^ in terms of the distance d{x) from Ig. Grieser and Jerison (1998) study 
the first eigenfunction {k = 1) in detail. The approximation (3.4) is shown to hold 
on and near I( with an error 0{s/i), and, if loc. max. / denotes the point where a 



it is not hard to show that the inradius is s + o(s) tlicii 

^This comphcates the analysis considerably since when applying A to (3.4) then one obtains, 
in addition to the terms yielding Pquo, terms involving first and second derivatives of h and hi, 
and the point is to control these. 

the cited papers only k = 1,2 are considered, but at least for the eigenvalues the results 
extend to arbitrary k easily. 

■'^"'^Note that different domains can have the same height function h. This question may 
be understood in two ways: First, one uses only 'elementary' (geometric) information about h. 
Second, one allows to use transcendental information, for example the eigenfunctions uq of the 
Schrodinger operator (3.3). We consider both questions. 
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function / assumes its maximum (this is unique for the first eigenfunctions ui , ipi 
in the case of convex domains) then 

(3.10) loc. max. = (loc. max. ipi,0) + 0{e). 

Grieser and Jerison (1996) study the second eigenfunction, in particular the nodal 
line, in detail. First, the approximation (3.4) is quantified to hold with an error 
0{e'^/i) near the nodal hue (eqn. (15) in loc.cit.), and this is used to show that the 
nodal line is contained in a vertical strip of width 0{s^). All these estimates are 
sharp in order of magnitude. See Subsection 3.3.1 for more on this. Jerison (2000) 
studies the first nodal line in the Neumann problem from a similar perspective. 

3.2.3. Curved and twisted tubes. Given a smooth embedded curve 7 : / ^ R" 
parametrized by arc length, a bounded open set Y C R"~^ and a family {^x)x£i 
of hnear isometries $3, : M."~^ ^-yix)! (the normal space to 7 in the point 7(0;)) 
depending smoothly on x we call the set 

f2e = {l{x) + e^^{y) : x&I, yeY} 

a thin curved tube with cross section Y. It is embedded for small e. We call it 

non-twisted if $ can be chosen such that for each y the curve x ^x{y) cuts each 
normal slice j{x) + Ran ^x orthogonally. It is easy to see that in this case the 
(Euclidean) Riemannian metric on O^, expressed in the coordinates {x, y) G I xY, 
is diagonal of the form^^ 

(3.11) g = di&g{a'^,e^,...,e'^), a{x,y) = 1 - e'y"{x) ■ ^x{y)- 

For example, any curved tube with Y a disk centered at zero is non-twisted; also, 
given 'generic' 7 and Y, there is a non-twisted curved tube, unique if fixed at one 
end. See Preitas and Krejcifik (2007) and the references given there. 

What should one expect for the eigenvalues of a curved, non-twisted tube with 
constant cross section, with Dirichlet boundary conditions? The Laplacian on fi^ 
expressed in the coordinates (x, y) is the variable coefficient operator^^ 

n-l 

Ag = -a~^dxa~^dx - e~^a~^ ^ ^y^ady^ 

3=1 

which is symmetric on its Dirichlet domain in L^{I x Y,adxdy). Then P = 
a-^/^Agfl"-^/^ is a unitarily equivalent operator in L^{IxY, dxdy). Symmetry implies 
that it has no first order terms, and indeed from (3.11) one sees that 

(3.12) P = -a-'^dl + e-^Ay + V, 

where Ay is the standard (Euclidean) Laplacian on Y and Ve{x,y) — — |k(.t)^ -t- 
0{e) is smooth in s,x,y. Here k{x) = \y{x)\ is the curvature. Since a = 1 + 0{s) 



l^Writc xo = X, (xi,. . .,Xn-i) = y. Then g^j = dx,'^' ■ dx^'S' where ^{a;, j/) = yix) +s^x(:y)- 
Since is an orthogonal Unear map for each x, we have gu = for i > 1 and gij = for i,j > 1, 
i^j- By construction 7' (a;) ± #2,(1/) and by the non-twisting assumption J- dy^^xiv) 

for all x,y. Therefore goj = (7' + edx^) ■ dy-^ = 0. Finally, by non- twisting, 7' and dx'^ are 
parallel, so goo = II7' = (1 + ■ dx^)^, and differentiating Yix) ■ ^x{'y) = in x yields 

'y'dx^ + 7"$ = and hence the given expression for goo = €? ■ 

""^■^Use the general formula A„ = — ,} dx^ g'-* Vdet gdx^ , where (g'-') is the inverse matrix 

V Get 9 ' J 

to g, and the volume element \/ det gdxo ■ ■ ■ dxn-i- 
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Figure 2. A rectangle with an end 

this is a 0(e) perturbation of — jI^'^)+£^^^y , which is a direct sum operator 

on / X y. Therefore, one expects the eigenvalues of P, and hence of A on il^, to be 

(3.13) Afc(e) = e^^t'o + /ifc + 0(e), /ifc = eigenvalues of - ^ —. 

This is precisely the result of Freitas and Krejcifik (2007) (see also the nice survey 
by Krejcifi'k (2008)). They also consider slightly more general operators of the form 
(3.12). In addition they show (assuming dY is smooth) that the eigenfunctions are 
0(e)-close to the product type eigenfunctions ^k{x)ipo{y / e) (with ipk the eigen- 
function corresponding to /i^), uniformly (!) with error linearly decaying at the 
boundary, as well as certain uniform derivative estimates. They use this to show 
that the nodal set (zero set) of the fc-th eigenfunction on il^ must lie within 0(e) 
of the nodal set of tp^ and that it must touch the boundary for small e. 

Curved and twisted tubes were studied in many papers, for example by Exner 
and Seba (1989a), Duclos and Exner (1995), Ekholm, Kovarik, and Krejcifik (2005), 
Bouchitte, Mascarenhas, and Trabucho (2007), mostly as waveguides. See also 
Geer and Keller (1983) for a (non-rigorous) WKB-analysis of curved waveguides 
and tubes with variable cross section. 

3.3. Intervals and curves as limits: Tubes with ends. What difference 
does an end make? Let us see what to expect, by adding an (ae) x e rectangle as 
an end to a 1 x e rectangle, for a fixed a > 0. This yields a (1 + ae) x e rectangle, 
so the eigenvalues change from c^^Tr^m^ + vr^P to e^^Tr^m^ + tt^P{1 + ae)^^ = 
e^^TT^m^ + vr^P — e 2a7r^P + O(e^), that is, the end makes a contribution of order 
e^, which is rather small and hence negligible for many problems. 

However, ends do make an important difference in a number of situations: 

• In some spectral geometry problems it is precisely the e^ term of the 
eigenvalue asymptotics which matters, so the ends matter. See 3.3.1. 

• In the global analysis context the influence of the ends is actually the 
main focus (see the explanations in Section 2). It is also stronger since 
the spectral invariants involve all eigenvalues simultaneously. 

• As we will see in 3.4.2, for graphs the precise form of the ends may influ- 
ence even the e~^ and e'^ term in the eigenvalue asymptotics. 

3.3.1. Rectangles with an end. Grieser and Jerison (2007) consider the Dirichlet 
problem on rectangular domains with one end added: Fix a Lipschitz function 
g : [0, 1] [0, oo) and set 

(3.14) n,^{{x,y)£R^;y^ey, < y < 1, -eg{y)<x<l}. 

See Figure 2. For the precise description of the eigenvalues and eigenfunctions, one 
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needs to consider the rescaling (1.5) and its limit as iV ^ oo, defined as 

(3.15) = {{x,y) e M^; < y < 1, -g{y) < x}. 

This is a half infinite strip with an 'end' added at the left side. The Laplacian 
on X°° has continuous spectrum [77^,00). We need some information^^ about the 
spectral value tt^: 

There is a unique function U on X°° satisfying AU = n'^U in 

C/ = at dX°° and U{x,y) = xsm{wy) + (9(1). 
Furthermore, there is a unique number a = a{g) > such that 

(3.16) U{x,y) = {x + a{g))s\n{T:y) + 0{e-^), x ^ 00. 

a{g) is related to the scattering phase, sec footnote 31. (3.16) may be interpreted 
as saying that, spectrally, X°° is similar to the half infinite strip (— a(g), 00) x (0, 1) 
with Dirichlet boundary conditions. Therefore, one may expect that the spectral 
data for X'^ should be similar to that for the rectangle {—a(g),N) x (0, 1). This 
is in fact true to rather high precision. Grieser and Jerison (2007) prove (scaling 
back to Oe) 

(3.17) A.(£) = e-V + ^^^^+0(£3) 

(3.18) Uk,s{x,y) = sin(fc7r^-^-^) sin(£:~^7ry) + 0{e^) for x > Seloge"^. 

for the suitably normalized fcth eigenfunction, where the estimates are uniform (for 
fixed k) with all d^, dy derivatives. 

Here is a question where a precision of estimates that 'detects' a is needed. 
Consider the problem whether the estimate (3.10) is sharp in order of magnitude, 
i.e. whether 0{e) may be replaced by 0{£^) for some r > 1. To prove a negative 
answer requires constructing an example which is far enough from product type 
but still can be analyzed sufficiently explicitly. This can be answered using (3.18), 
and it is the quantity a{g) which is essential here: (3.18) implies that loc. max. ui = 
\ — + 0(£^). So if one finds two functions 51,(72 inducing the same height 
function h but having a{gi) ^ 0(52) then it follows that (3.10) is sharp. This is 
carried out in (Grieser and Jerison, 2007). 

3.3.2. Global analysis. Since the emphasis in this article is on methods common 
to the fields mentioned in the title, and in order to avoid having to introduce a lot 
of notation, we only mention one recent result from global analysis. 

We need to recall the definition of the determinant of a self-adjoint, non- 
negative elliptic operator P, for example the Laplacian. Let Ai < A2 < . . . be 
the positive eigenvalues of P. Set C,p{s) = '^JLi ^J'^- This converges and is holo- 
morphic for s G C with sufficiently large real part, and can be shown to have a 
meromorphic continuation to s € C which is regular at s = 0. Then define 

detP:= e-f^(°^ 

This is motivated by the fact that the analogous formula is true for invertible linear 
operators on a finite-dimensional space, with the usual determinant. 



^■^It is also important that there are no eigenvalues < vr^ (sec the explanations in 3.4). 
This is true (and easy to see, using the variational principle) in this setting, because the end is 
contained in the doubly infinite strip M X (0, 1). 



14 



DANIEL GRIESER 



Consider a tube of the form (1.4) with / = [—1, 1] and assume it has smootli 
boundary. We state the result in terms of the corresponding long cylinder 
defined by (1.5). This consists of the ends X(, Xr and a cylinder [— A^, A^] x F con- 
necting them. Miiller and Miiller (2006) analyze (among other things) the asymp- 
totics of detP^Jv as N ^ oo, where P^n is a Laplace type operator of product 
type. Laplace type means that it is a self-adjoint differential operator whoso loading 
part is the Laplacian, and product type means that it equals — + Py over the 
cylindrical part [—A/', A?'] x Y, for a Laplace type operator Py on Y}^ If X^ has 
boundary then Dirichlet boundary conditions are imposed. 

Let X^ = [Xi U ([0,oo) X Y)] /ji, i.e. the left end with a half infinite cylinder 
attached at the boundary Y, and define X^ similarly. These are the two compo- 
nents of the limit space introduced below in (3.21). Since Pjc" is product type over 
the cylinder part, it induces operators Px^ , Px^ on Xf° , X^. Assuming that Py, 
Pxf° and Px^ are injective (for example, P = A -|- 1), Miiller and Miiller (2006, 
Theorem 1.2) prove 

det PxN ~ e^^^P, N ^oo 

for certain constants A, B. A is calculated in terms of the function (^p^ , and B in 
terms of similar quantities (relative determinants) on and det Py. 

There has been a lot of work in similar directions, sec the references in Afiiller 
and Miiller (2006). As examples we mention Park and Wojciechowski (2005), Loya 
and Park (2004), Lee (2006). An analysis of analytic torsion, involving similar 
issues, was done by Hassell (1998). 

The basic principle and spectral invariants: det P, like most spectral invariants 
occur ing in global analysis, depends on all eigenvalues of P simultaneously. There- 
fore one might expect that the basic principle, which deals with Afe(£) for fixed 
k, is not applicable. However, in singular limit problems one can often separate 
the influence of 'small' eigenvalues from the influence of 'large' eigenvalues. In the 
analysis of the former the basic principle does apply, and therefore some of the 
methods in Section 4 are common to fixed k and spectral invariant analysis. For 
the large eigenvalues one uses other techniques, for example heat kernel asymptotics 
as f — > 0. 

3.4. Graphs as limits. In this section we consider tubes modeled on a finite 

metric graph G = flo- We use the notation introduced in Section 1. We first 
introduce quantum graphs and then consider the following situations: 

• General boundary conditions, constant cross section along all edges; 

• Neumann boundary condition, cross section may vary; 



"'^'^Morc generally, is allowed to be a matrix- valued operator (or system of operators); 

that is, a vector bundle E over is given, and P^jv acts on the sections of E. E and P^a- are 
assumed to be of 'product type' over the cylindrical part of X'^; that is, E is the pull-back of a 

,2 

bundle Ey over Y and Pxn = + for a Laplace type operator Py acting on sections of 

Ey over Y. Laplace type means that, with respect to a (product type) hermitean structure on 
E, the leading part of Pxn, when representing it in terms of an orthonormal frame of E, is Aj,^ jv 
times the identity on E. 

l^That is, the situation described in Section 1 with all Ye the same. The case of different 
Ye is also treated by Grieser (2007). The main difference is that edges e for which the lowest 
eigenvalue uo(Ye) is bigger than min^/g^ ^o{Ye') should be deleted in the quantum graph of the 
main result in 3.4.2. 
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Thin tubes modeled on graphs have been studied extensively in the last decade. 
The case of Neumann boundary conditions (or empty boundary) has been well- 
understood for some time. For Dirichlet or more general boundary conditions, 
decisive progress was made only recently. We will explain some of these works. 
First, we introduce the limit operators on the limit graph. 

3.4.1. Quantum graphs. What should one expect for the eigenvalues and eigen- 
functions on a tube around a graph? 

In the two situation mentioned above, the first cross-sectional eigenvalue is 
constant (say vq) along each edge {vq = in the Neumann case, with the constant 
function as eigenfunction; Pq possibly positive but constant in general). From the 
basic principle, see Section 3.1, one expects the fcth eigenfunction Uk^s on to be, 
in the tube around any edge e, close to a product u'f,{x)ipQ{y /e), where u'f, is an 
eigenfunction of on the edge. Also, the eigenvalue should be approximately 

e~'^uo + /Ufc, where /Xfc is the eigenvalue corresponding to Now it must be the 
same /ife for all edges e, since Uk^e is an eigenfunction on all of Qg. If one knew how 
the tt^ fit together at the vertices then we could determine the ^k- Without such 
information the cannot be determined. This 'fitting together' is made precise in 
the concept of a quantum graph. 

Thus, a quantum graph is, per definition, a metric graph together with a self- 

.2 

adjoint realization of the operator — ^ along each edge; such a realization is given 
by boundary conditions at the vertices. One may characterize all possible such 
boundary conditions, see Kostrykin and Schrader (1999). We will only encounter 
a special class of boundary conditions described as follows: 

For each vertex v denote the set of adjacent edges by E{v) = {e € E : e ^ w}. 
Suppose for each v a linear subspace 

is given. This determines a boundary condition as follows: Suppose a function 
is given on each edge e (including its endpoints). We say u = {u'')eeE satisfies 
the boundary condition given by the iff for all v 

(3.19) (u"(«))e~. e 

(3.20) (a„?i^(«))e^. e w^^. 

Here, u"^ is the restriction of u to the edge e, dnu'^{v) is the derivative of u'^ at v 
pointing inside e, and W:^ is the orthogonal complement of Wy. Two special cases 
are especially important: 

Kirchhoff boundary conditions: All Wy = span{(l, . . . , 1)}, that is, u is contin- 
uous at the vertices and the sum of derivatives at each vertex is zero. 

Decoupled Dirichlet conditions: All Wy = 0, that is, is zero at both ends of 
e, for each edge e. 

For more information on quantum graphs, see the recent surveys (Kuchment, 
2004) and (Gnutzmann and Smilansky, 2006). 

3.4.2. General boundary conditions, constant cross section. We assume that 
boundary conditions on dQ,e are arbitrary except that they respect the product 

decomposition [0, Z(e)] x eY^ along each edge. That is, the boundary of each cross 
section Yg is divided into two parts 2?, A/" and the parts of the vertex neighborhoods 



One may also consider more general quantum graphs, for example allowing a magnetic or 
an electric field. 



16 



DANIEL GRIESER 




Figure 3. Rescaled graph neighborhoods: X and X 



Xy where no edges are attached are decomposed similarly. This induces an anal- 
ogous decomposition of dn,;. Then we impose Dirichlet conditions on the V part 
and Neumann conditions on the Af part.^^ 

It turns out that the expectation stemming from the basic principle and ap- 
plied in the previous section is almost correct, and that the quantum graph (i.e. 
the boundary conditions at the vertices) can be determined in terms of certain 
scattering data. To state this, we define the limit X°° of the rescaled tube (1.5) 



(3.21) X°° ■.= [jX^ (disjoint union) 

V 

where the 'star' X^ of a vertex w of G is obtained by attaching a half infinite 
cylinder [0, cxd) x Ye to the vertex manifold Xi,, for each edge e incident to v, 
see Figure 3. Fix a vertex v. The Laplacian —Ax^ (with boundary conditions 
stemming from those on fl^) has continuous spectrum [^"0,00), and possibly L^- 
eigenvalues. Standard scattering theory produces for each A G [I'coo) a unitary 
scattering matrix 

(3.22) S„(A) : C^(") ^ C^(") 

describing the generalized eigenfunctions for the spectral value A, i.e. the solutions 
u of 

(3.23) Au = Au, 3N \u{x,y)\ = O(x^) for x-^ 00. 

For the determination of the quantum graph only the scattering matrix at A = i/q 
is needed. Using separation of variables it is easy to see that any solution of (3.23) 
with A = 1^0 is, along each edge cylinder [0, cxd) x Y^, asymptotically equal to either 
''^efoiy) or WeXipo{y) for some We € C, where x is the coordinate along the cylinder 
axis. Then Sy{i>o) is characterized by the property that it is a unitary involution, 



"Robin boundary conditions may also be considered, see (Grieser, 2007). 
In the case of (3.14) this would be the union of (3.15) and a half strip [0, 1] X [0, oo), so 
the notation is not consistent with (3.15). However, the two ways of looking at the 'limit' of 
are equivalent, see also footnote 30. 
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i.e. it has only eigenvalues ±1 with mutually orthogonal eigenspaces, and that its 
+l-eigenspace is given by the bounded solutions, i.e. 

(3.24) Sv{i'o)w = w there is a generalized eigenfunction u{x,y) wipo{y) 

Here, ii w = {we)e<£E{v) then wipo is the function that equals We^Po on the edge e. 
Main result on spectral convergence for tubes modeled on graphs: Let 

Hk be the eigenvalues, counted with multiplicity, of the quantum graph given by 
the metric graph G with boundai-y conditions (3.19), (3.20) where Wv is the +1- 
eigenspace of Svii'o), for each vertex v. Also, let n < T2 < ■ ■ ■ < td be the 
L'^ -eigenvalues < of /^x^^- Then the eigenvalues on Oe satisfy, for fixed k and 
e — > 0, and for some constant c > 0, 



This result has features distinct from the special case of Neumann boundary 
conditions, which had been known before and which will be discussed below: 

• The eigenvalues (3.25) lie below those expected from the 'basic princi- 
ple'. Their presence (i.e. _D > 0) is fairly typical in the case of Dirichlet 
boundary conditions, as is well-known in waveguide theory, see for exam- 
ple (Schult et al., 1989), (Avishai et al., 1991) and (Duclos and Exner, 
1995). 

• The boundary conditions at the vertices of G are determined by scattering 
data. 

We will see in 3.4.3 that the behavior of the Afe(e) in the Neumann problem is 
determined by the combinatorics of G alone. These features show that this is not 
the case for the Dirichlet problem; here the particular geometry of the edge and 
vertex neighborhoods matter also. 

This result is proved by Grieser (2007), and partially by Molchanov and Vain- 
berg (2007), by a different method (see Section 4). Molchanov and Vainberg (2007) 
also admit graphs with infinitely long edges (see also Molchanov and Vainberg 
(2006a), and Molchanov and Vainberg (2006b) for a simpler case). Grieser (2007) 
derives a full asymptotic development of Afc for all k (cf. (3.17) and the comments 
before and after it for reasons to consider at least a third term) , even uniformly for 
k < Ce~^, for some constant C,^° and obtains asymptotics for the eigenfunctions 
as well. If the vertex neighborhoods are 'small' in the sense that their Dirichlet 
eigenvalues (however, with Neumann condition at the parts of the boundary where 
the edges are attached) are bigger than i/q then D = and all Wy = 0. In this case, 
the theorem was proved by Post (2005). 

A different m,odel: Although the emphasis is on bounded tubes in this article, wc 
mention interesting recent work on certain unbounded tubes modeled on the graph 
SIq having two infinitely long edges joined at one vertex, with Dirichlet boundary 
conditions: Given a non- intersecting smooth curve 7 : R ^ with compactly 
supported curvature k (so it consists of two straight half lines outside a compact 
set) define fie as the £-neighborhood of the curve 1 1— > s^'y{t), for some r < 1. Thus, 
the curved portion (corresponding to the vertex neighborhoods in the graph model) 



(3.25) 
(3.26) 



Afe(£)=£-Vfc + 0(e-^/^), 
Afe(e) = e~'^i'o + Hk-D + 0{s) 



k = l,...,D 



k>D 



this corresponds to analyzing the full interval [1^0, ^^i) ™ the long cylinder problem 
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has length ~ e^, and the same thickness e as the rest of the tube. Albeverio, Cac- 
ciapuoti, and Finco (2007) prove, assuming r < 2/5, convergence of the resolvent 
(in a certain sense) to the resolvent of the operator —dP/dx^ on ^q, with one of 
two possible boundary conditions at the vertex (either Dirichlet decoupled or with 
a certain coupling). Which boundary condition occurs depends on whether the 
operator —-j^ ^ 1^"^ ^ bounded generalized eigenfunction at the bottom (zero) 
of its continuous spectrum. This is very similar to the situation in the main result 
above: The limit operator is the graph Laplacian, with boundary condition at the 
vertex determined by the rcscaling of fi^ which makes the vertex neighborhood of 
a fixed 'size'. Cacciapuoti and Exner (2007) discuss a slightly more general model 
(with the same SIq) in which they can produce more general boundary conditions 
at the vertex. 

3.4.3. Neumann boundary conditions. In the case of Neumann boundary con- 
ditions on all of dflg it is easy to see that, in the main result of the previous section, 

D = and the quantum graph boundary conditions are Kirchhoff at each vertex. 
Also, fo = 0, and therefore the 'interesting' behavior occurs already in the lead- 
ing term. The result Xk{s) — > /Xfc (e — > 0) is then accessible by more elementary 
methods (quadratic forms) and is also much more stable than the result for gen- 
eral boundary conditions, for example under bending of the edges or varying the 
thickness of the tube along the edges. Such results^^ were obtained much earlier. 

Colin de Verdierc (1986) seems to be the first who proved this result (for straight 
tubes of constant width) , but it seems that the paper was not known in much of 
the mathematical physics community until recently. 

Freidlin and WcntzcU (1993) take up the problem from a probabilistic per- 
spective; they do not consider the spectral problem but analyze the equation 
(A + c)u = g on. (a neighborhood in of CIq) for given continuous func- 
tions c > and g and show that the solution u converges for £ ^ to the solution 
of the corresponding equation on fio, with Kirchhoff boundary conditions at the 
vertices. This is, in essence, strong resolvent convergence. Saito (2000) obtains a 
similar rc^sult for curved tubes in the plane of variable thickness. 

Rubinstein and Scliatzman (2001) prove the spectral convergence for curved 
tubes in the plane of variable thickness. They also allow a magnetic field. Their 
method was simplified and the result gc^neralized to quite arbitrary vertex neigh- 
borhoods (of 'roughly the size' e) by Kucliincnt and Zcng (2001), and extended to 
vertex neighborhoods that shrink more slowly (at a rate e'', r < 1) by Kuchment 
and Zeng (2003). Depending on r this may give different boundary conditions at 
the vertices. Exner and Post (2005, 2007) prove similar results for straight tubes 
of constant cross section in the manifold setting. 

4. Methods 

What makes thin tube problems difficult? There are two sources of difficulty 
when analyzing spectral problems of the Laplacian on fl^, uniformly as £ — > 0: 

• In a tube without ends, one has to deal with two scales: The transversal 
direction scales like e while the longitudinal direction is fixed. This is 
closely related to so-called semi-classical problems, i.e. the analysis of the 
operator (3.3) as £ — > 0. 



often with weaker error terms than 0{e) 
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• In a tube with ends one has the additional difficulty that different parts 
of ris scale in different ways: the ends scale like e in all directions, the 

tubes around the edges only in the transversal direction. 

There are various ways to deal with these problems. They differ in complexity, 
flexibility and precision. We mention some of them. They will be explained in more 
detail below. 

• Quadratic forms: This is the simplest and most flexible tool. It is quite 
robust in low regularity situations. For example, it allows a fairly direct 
and uniform treatment of tubes around graphs where the vertex neigh- 
borhoods scale at powers with r ^ 1, or where they are not given as 
scalings of a fixed model. 

However, the quadratic form technique does not give optimal results 
in many situations, especially when one needs to understand lower order 
asymptotic terms. For example, it has not been successful for the Dirich- 
let problem on tubes with ends (or, more generally, tubes modeled on 
graphs). 

• Matched asymptotic expansions: This is a classical tool for problems in- 
volving regions that scale in different ways. It can give very precise infor- 
mation. A technical problem is to show that the solutions (eigenfunctions) 
constructed by this method are all the solutions. 

• Resolvent analysis: A good understanding of the resolvent of an opera- 
tor yic;lds information on its spectral data and also much more. Rough 
information on the resolvent may sometimes be obtained using quadratic 
forms. A precise description of the resolvent is rather complicated since 
in addition to the parameter s one has to analyze the dependence on the 
resolvent parameter A. There are highly sophisticated techniques to do 
this, and when applicable they yield the most precise information. The 
pseudodifferential technique described below allows to treat tubes whose 
cross section varies mildly. 

• Matching of scattering solutions: This is a technique especially adapted 
to thin tube problems with constant cross section. It is rather direct and 
yields very precise information. As for matched asymptotics, one has to 
work a little harder to show that one obtains all eigenfunctions with this 
method. 

4.1. Quadratic forms. The simplest method to obtain eigenvalue compari- 
son results is to compare the quadratic forms corresponding to the operators in- 
volved: If P is a self-adjoint operator with domain Dom(P) and quadratic form 
qp{u) = {u, Pu) then the fcth eigenvalue of P is (assuming P has discrete spectrum 



■^^ Another, classical example for this weakness of quadratic forms occurs with the Weyl 
formula for the asymptotic distribution of eigenvalues of the Laplacian on a bounded domain 
Q C K" with boundary: Let JV(A) denote the number of eigenvalues < A. With Dirichlet- 
Neumann bracketing, a quadratic forms technique, one can prove 7V(A) = C„A"/2 -I- R(A) with 
an error estimate -R(A) = 0{A("^^-'''^ log A), see Courant and Hilbert (1924), for example. The 
sharper estimate -R(A) = 0(A'"^-'^'/'^) holds but apparently cannot be proven using quadratic 
forms. There are various proofs of it using resolvent or heat kernel analysis. See Grieser (2002) for 
further references. On the other hand, the Dirichlet-Neumann bracketing still works for domains 
with very rough, for example fractal, boundary, see (Lapidus, 1991). 
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and is bounded from below), by the variational principle, 
(4.1) Afc(P) = inf _ sup 

where Lfc runs over fc-dimcnsional subspaces. This immediately implies that if 

Qp' < Qp then \k{P') < ^k{P) for all fc, and more generally: 

Comparison principle: If J : Dom(P) Dom(P') is a linear 
map that decreases the norm at most very little and increases 
the quadratic form at most very little (i.e. qp'{Ju) < qp{u)+ 
small) then Xk{P') can be at most slightly bigger than Afc(P). 

Here, 'very little' must be measured in an appropriate norm (e.g. the Sobolev 
norm if P, P' arc second order differential operators 



s 
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This principle allows for P, P' to live in different Hilbert spaces and is therefore 
suitable for our problem, since we want to compare A^, defined in L'^{Q.s), with 
Ao, defined in L'^{Q.q). The quadratic forms in our context are simply qA^{u) = 
In. |V«P. 

We explain how the comparison principle yields the eigenvalue convergence in 
the case of thin tubes around a graph, with Neumann boundary conditions, see 
3.4.3 (but with a worse error 0{^/e) in (3.26); we loosely follow (Exner and Post, 
2007) which we refer to for more details, and to (Exner and Post, 2005) for the 
full proof along these linos). We need to find operators J : L'^ifla) L'^{fle) and 
J' : L'^{ile) i^(f2o) respecting the domains of A which yield good comparison 
estimates of norms and quadratic forms as needed in the comparison principle. 
Good estimates may be expected if J, J' reflect the expectation that eigenfunctions 
on rig with fixed k and e small should be approximately constant in the transversal 
direction. 

To produce from a function / on the graph Qq a function J(/) on is easy: 
Simply extend / to be constant transversal to each edge and constant in each vertex 
neighborhood, then multiply by e""*/^ vol(F)~^/^ to make up for the area of the 
transversal slice. It is then easy to obtain the needed estimates, and this yields 
Afc(A,)<Afc(Ao) + 0(£). 

The opposite direction is a little trickier: Given a function u on fie, it is natural 
to define the function J'{u) on flo along the edges by taking transversal averages, 
and at each vertex by taking an average of u over the vertex neighborhood (and 
then multiplying by a suitable scaling factor). However, the resulting function may 
be discontinuous at the vertices, so some kind of 'fudging' has to be done. For this, 
various procedures have been proposed in the literature, the simplest one being the 
use of a smooth transition function^^. Then the estimate qo{J'{u)) < ge(u)+'smair 
is fairly straightforward. The hardest part, however, is to show that ||'/'(w)||x,2(Q(,) 
cannot be much smaller than ||u||i2(Q^'). This is done in two steps: 

(1) Consider a transversal slice to a point of an edge. Then the deviation (in 
of the slice, times the scaling factor) of u from its average over the 



■^^See, for example, Exner and Post (2007) for a precise statement, or Colin de Verdiere 
(1986) for a different version. 

■^■^For example, if /o is a given continuous function on [0,/] and a £ M is close to /()(0) then 
one may produce a function on [0, 1] which is equal to a at zero and close to /o everywhere by 
choosing a smooth p : [0,1] — > [0, 1] which equals one at a; = and zero for x > 1/2, and setting 
f{x) = fo{x) + p{x){a-fom. 
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slice is 0(e)||'u||_f/i(siice) by the Poincare inequality^^. Integrating along 

the edges yields a lower bound of || J'(w)||L2(no) by the L^-norm of u over 

the edge tubes (minus a small error 0(e)||?ij|^i(n_)). 
(2) Bound the L'^-norm of u over the vertex neighborhoods by a small multiple 

of its i^-norm over the edge tubes (plus small ||u||iji -contributions). 
Together, this yields the bound |1 J'(m)1|l2(s2„) > \\u\\l2(^q^)-0{^)\\u\\hi{q,), hence 
Afc(Ao) < \k{A^) + O(y^), which was to be shown. 

It is in Step 2 that the proof hinges most essentially on the fact that Neumann 
boundary conditions arc considered. It can be still salvaged for Dirichlct conditions 
provided the vertex neighborhoods are 'small', as observed by Post (2005). In 
general however, a bound as needed here cannot hold, due to the existence of L^- 
solutions to the associated scattering problem, see footnote 32. In fact, will then 
have eigenvalues that do not come from any graph eigenvalues (those in (3.25)). 

The quadratic forms technique was also used by Colin de Verdiere (1986), 
Rubinstein and Schatzman (2001) and Kuchment and Zeng (2001, 2003). 

4.2. Resolvents, I. The resolvent Rp{z) = (P — z)~^ of an operator P, de- 
fined for z e p{P) = C \ cr(P), is a very powerful object in the study of spectral 
degeneration. If one can describe it sufficiently well, then one may deduce infor- 
mation about the eigenvalues and also the eigenvectors: if P is self-adjoint with 
discrete spectrum and 7 is a simple closed curve in p{P) then 

(4.2) trace / zRp{z)dz, / Rp{z)dz 

are the sum of the eigenvalues of P inside 7 and the orthogonal projection to the 

corresponding spectral subspacc, respectively. Therefore, if we have a family P^ 
and 7 C p{Pe) for each e then the £-dependence for these objects can be deduced 
easily from that of Rp^. In the case of regular perturbations this was studied in 
detail by RcUich (1936), see also Kato (1976). 

Denote by R^ the resolvent of the Laplacian on fi^ (in any one of the situations 
discussed above) and Rq the resolvent of the (expected) 'limit' operator on Oq- Since 
degeneration of a tube to a one- dimensional limit object is a singular perturbation, 
the behavior of Rp^ as e ^ is more complicated than for a regular perturbation. 

There are two approaches to understanding R^: 

(1) Compare R^ and Ro directly, or 

(2) construct an approximation (parametrix) for Rg,^^ then show that it is 
actually close to R^. 

Direct resolvent comparison: In this approach, an initial difficulty is that Re 
and Rq act on different Hilbert spaces. In the case of a tube without ends (see 3.2.1 
for notation), one may write L^{fle) = Lq(B L>i where Lq corresponds to the first 
transversal mode and L>i to the higher ones. Lq may be identified with L^{I) = 
L^(rJo), and therefore i?o niay be regarded as operator on Lq, and extending it 
by zero on L>i one obtains an operator on L^{^lg) which may be compared with 

This is the approach taken by Priedlander and Solomyak (2007a,b), except that 



25The Poincare inequality says that /y \u - avey up < Cy fy |Vup where avey u is the 
average of u over Y and Cy is a constant depending on Y (but not on u) . A change of variable 
gives the scaling law Cgy = s^Cy ■ 

■^^this may, and probably will, involve some information coming from Ro 
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they subtract the divergent term e^^tt^ from Aq^ and (3.3) before taking inverses. ^'^ 
The closeness of these operators is proved by a quadratic form comparison. 

In the case of tubes with ends, or more generally graphs, a similar program 
was carried out by Post (2006) (see also Exner and Post (2007, section 4)) in the 
case of Neumann boundary conditions. However, so far this approach has been 
unsuccessful when dealing with the Dirichlet problem for tubes with ends, or more 
generally tubes modeled on graphs. 

Resolvent construction: Molchanov and Vainberg (2006a, 2007) use scattering 
solutions (sec (4.8) below) to construct a candidate approximation for the integral 
kernel of for a tube modeled on a graph, with Dirichlet boundary conditions. 
They outline an argument implying that this approximates Rg to order 0{s) in 
operator norm. A construction of R^ to all orders was given by Hassell, Mazzeo, 
and Melrose (1995). The method of the latter will be explained in 4.4. 

4.3. Matched asymptotic expansions. This is a general technique for prob- 
lems involving different regions (sometimes called regimes) that scale in different 
ways. It was used by Grieser and Jerison (2007) in the setup introduced in 3.3.1. 

The basic steps are: 

Step 1: We would like to find an eigenvalue A(£) with eigenfunction Ue on fig. One 
assumes that the eigenvalue has a complete asymptotic expansion as £ — > 

oo 

(4.3) A(e) ~ e-^TT^ + ^ eU^*) , A^'^ G R. 

i=0 

For the eigenfunction one assumes expansions of two kinds, corresponding to the two 
regimes of ilg. The first regime is the interior of the edge: For any fixed x G (0, 1] 
and y G [0, 1] one assumes an expansion refining (3.4) (with h= 1, hi = 0) 

oo 

(4.4) Ue{x,ey) ^^e'u'''\x,y) 

i=0 

for functions m^*^ on ^lQ := (0, 1) x [0, 1]. The second regime is the order e neigh- 
borhood of the left end: For any fixed {x,y) G X°° one assumes an expansion 

oo 

(4.5) Ue{£x, ey) s-'v^-'^ {x, y) 

for functions v^^"^ on 

Step 2: Plugging first (4.3), (4.4), then (4.3), (4.5) into the equation Awg = X{e)ue 
and comparing coefficients of powers of e one obtains two infinite sequences of 
PDEs involving the coefficients A('\u/'',f(^'. Since satisfies Dirichlet boundary 
conditions at cJfie, the u^^^ vanish at all sides of f^o except possibly at a; = 0, and 
the t;^-'^ vanish at dX°°. This is not enough information to solve the equations 
uniquely. In addition, one needs matching conditions coming from the fact that the 
large x behavior of the v^^^ should be related to the behavior as a; ^ of the u'^^\ 
These conditions are obtained by setting x = ex in (4.4), expanding each term in 
Taylor series around a; = and comparing coefficients of e powers with (4.5). 



^^Here, only inverses are needed, not the values of the resolvent at other spectral parameters. 
The closeness of eigenvalues is deduced from the simple estimate \iJ.j(A) — ij,j{B)\ < ||j4 — _B|| for 
compact selfadjoint operators, applied to these inverses. This is simpler than (4.2). 
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The PDEs, with these boundary and matching conditions, can now be solved 
recursively. For each A*^*^^ — fc^Tr^, fc S N, they determine A^'^ for i > 1 and u'*-*, 
w^^^ for j > uniquely, up to a normalization factor for u^. 

Step 3: One can now construct an approximate solution by 'gluing' the formal 
expansions (4.4), (4.5) in a suitable way. Cutting off the series after the Mth term 
one obtains Us,k,M, ^k,M{s) satisfying the boundary conditions and 

(4.6) ^Us,k,M = A/s,M(e)w£,fe,M + O(e^). 

One needs to show that this implies that u^,k,M , ^k.Aiis) are close (to order e^) 
to an actual eigenfunction and eigenvalue. Also, one needs to show that all eigen- 
functions are obtained in this way. Both are quite straightforward in the setting of 
domains (3.14) since by domain comparison one has good a priori estimates of the 
eigenvalues of fig. Closeness of approximate to actual eigenfunctions follows from 
a spectral gap argument. This works for all M, so it follows that Afe(£) and Uk,s 
actually do have expansions as in (4.3), (4.4), (4.5). 

4.4. Resolvents, II: PseudodifFerential operators. The classical pseudo- 
differential operator (\1/D0) calculus provides a systematic way to study the re- 
solvent of an elliptic operator on a compact manifold X (see (Shubin, 1987) for 
example). More precisely, if R{X,p,p') is the integral kernel of the resolvent, where 
p,p' G X and A is the spectral parameter, then the standard \1/D0 parametrix 
construction gives an explicit way to find the asymptotic expansion of R as p ^ p' 
('at the diagonal'), to all orders. 

In the last two decades an intensive effort has been made (still ongoing) to ex- 
tend the ^DO calculus to settings involving either singularities (or non-compactness) 
of the underlying manifold or singular perturbation problems, sec for example (Mel- 
rose, 1991, 1996) or (Schulze, 1991, 2001). Here, the goal is to obtain the asymp- 
totics of R not only at the diagonal, but also as one or both points p, p' approach 
the singularities (or 'infinity'), or as the perturbation parameter approaches a sin- 
gular limit. R. Melrose and his school have developed a very geometric and general 
approach to this problem and applied it to many kinds of singularities and degener- 
ations. Sec (Grieser, 2001) for an introduction to these ideas and more references. 

Hassell, Mazzeo, and Melrose (1995) studied from this perspective the degenera- 
tion of a thin tube in the case of closed manifolds, with an interval as limit space. We 
cannot explain this here in detail but only give a rough idea. The resolvent kernel is 
a function R{e, X,p,p') on the space X = {{e, X,p,p') : e S (0, 1], A ^ M, p,p' G fig}, 
and the goal is to obtain full information of the asymptotics at the 'boundary' and 
at the diagonal of this space, i.e. as e ^ 0, A ^ (in this case, i^o = 0), p ^ p' 
(allowing any combination of some of these limits), uniformly over all p,p'. This 
involves a large number of different regimes, similar to but more complex than the 
situation in 4.3 where we considered different regimes for the eigenfunction as func- 
tion of e and p. The goal is accomplished in essence by a combination of the classical 
^DO technique with the method of matched asymptotic expansions. A difficulty 
lies in taming the high combinatorial complexity of all the regimes and resulting 
matchings. This is made feasible by encoding these geometrically in a manifold 
with corners X which is a compactification of X, whose minimal boundary faces 
correspond to the regimes and boimdary hypersurfaces to matchings, cf. (Grieser, 
2001, section 2.5.3). See (Hassell et al., 1995, Figure 2) for a nice picture of X. 
More precisely, the work proceeds by first defining a suitable space X (motivated. 
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among other things, by model calculations), then defining a class of operators (a 
'A'-pseudodifferential calculus') as those whose integral kernels arc distributions on 
X having certain 'nice' (polyhomogeneous conormal) behavior at the boundary of 
X and conormal (\E'DO) singularities at the diagonal, then showing that this class is 
closed under composition and has various other 'good' properties, then constructing 
a first order parametrix for A — A (this amounts to solving certain limit or model 
problems, one of which being the inversion of the Laplacian on an interval, another 
the scattering problem on X°°) and then improving this to any order by the stan- 
dard iteration in parametrix constructions. This can then be used to obtain full 
asymptotics of eigenvalues and eigenfunctions by (4.2). 

4.5. Matching of scattering solutions. This very effective method was 
used by Grieser (2007) for the proof of the main result in 3.4.2. We explain it 

in the simple case of a single edge graph with one end attached, as in (3.14).^^ 

We work on the rescaled domain X'^ = {(x,y) : < y < 1, —g{y) <x< N}, 
and its limit X°° defined in (3.15). We first need some facts from scattering theory 
on X°°, extending (3.16). For each a G (0, y/Sw) there is a unique bounded solution 

Eo, of 

(4.7) AEa = XEa on A = tt^ + a^, E^iox-- = 0, 
which has the form 

(4.8) E^ {x, y) = (e-^"^ + S{a)^°'^) sin ny + (x, y) 

for X > 0, where S{a) G C and Ra{x,-) is in the span of the higher transversal 
modes sinTrmy, m > 2, and consequently is exponentially decaying as x x.^^ 
S{a) is called the scattering matrix.'^" It is unitary and extends to a holomorphic 
function of a in a neighborhood of in C, satisfying 5(0)^ = 1. In our example 
S{0) = — 1 (this means that the +1 eigenspace in (3.24) is zero). Write 

(4.9) S{a) = e'f''"^ p holomorphic in a, p(0) = tt. 

p{a) is real for real a. It is called the scattering phase. 
Fix iV > 0. Suppose a is such that 

(4.10) e"*"^ + 5(a)e*"^ = 0. 

Then Ea{N,y) = Ra{N,y) = 0{e~'^^). Therefore, satisfies Dirichlet boundary 
conditions at dX^ , up to an exponentially small error. Also, AEa = XEa in X^ . 



■^^but the special structure of the end in (3.14) is inessential here, in particular need not 
be contained in the strip < y < e 

■^^This decay follows from Ra{x,y) = X^^=2 ''"'(s:) sin Trmy with rm{x) = 
2 /q Ra(x,y) sin nrriy dy, hence bounded, and satisfying — \/w?'k'^ — \Tm = from (4.7), 
which has only exponential solutions for A < 47r^. 

•"^ In this simple example it is a number, i.e. a 1 X 1 matrix. Note that in the general setup 
of (Grieser, 2007) this example would be treated as the union of two half strips, one emanating 
from each end (vertex), see (3.21) and footnote 19, and the scattering matrix would be the 2x2 
matrix diag(5(a), —1); it is easy to see from this that the two treatments are equivalent. 

''"'^The relation of (4.8) with (3.16), which is a scattering solution for a = 0, is this; Write 
p(a) = TT + aa{a) with <t(0) = p'(0). Then Ea{x,y) = -2ie*"''("'/2 gj^ j'^^,^ ^^)a] siuTry + 
Ra{x,y). Therefore the first transversal mode of a^^Ea converges, as a — » 0, to —2i(x + 
^^)sin7r5. It is easy to see that a~^i?a also converges. Therefore U = ^ lima_>o Q:~'^-Ea 
and a{g) = p'(0)/2. 



REFERENCES 



25 



From this it is easy to conclude that A hes within 0{e '^^) of a Dirichlet eigenvalue 
of A on 

Therefore, we need to analyze the solutions a of (4.10). Use (4.9) and introduce 
the variables z = aN, e = N~^, then this is equivalent to 

(4.11) p(ez) + 2z = 27rfc + TT for some fc S Z. 

By the implicit function theorem this can be solved locally for z, that is, for each k 
there is a unique analytic function s i— > Zk{s), defined near zero, so that z = Zk{e) 
satisfies (4.11) and Zfc(O) = kn. Then afc(e) = ezk{e) solves (4.10) with N = e^^, 
so one gets an eigenvalue TT^ +£2zfe(e)2 +0(e-^/^) on and hence an eigenvalue 

oo 

(4.12) s-\^ + + J2 4'^^' + 0(e-^/^) 

i=l 

for flg.^^ Also, -Eafc(e) is an approximation, with error 0{e~'^^'^), for the eigenfunc- 
tion. 

As in 4.3 one needs to show that one has obtained all eigenvalues in this way. 
For general graphs, this requires a rather involved argument, which essentially rests 
on a priori estimates for the eigenfunctions, see (Grieser, 2007). There are various 
pitfalls complicating the matter, for example the issue of multiplicities and almost 
multiplicities (eigenvalues can be very close together, so the spectral gap argument 
alluded to in 4.3 does not work). 

The method of matching of scattering solutions was also used in the global 
analysis context by Cappell et al. (1996), Miiller (1994), Park and Wojciechowski 
(2006). 
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